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ABSTRACT 


The  consequence  of  propellant  gas  flow  through  a  split  in  e 
cartridge  case  sidewall  is  usually  more  serious  with  aluminun  i;?.3es 
than  it  is  with  brass  cases.  In  order  to  gain  a  better  understa’...  ing  of 
the  failure  phenomenon  experienced  with  aluminum  cartridge  ies,  a 
combination  experimental  and  theoretical  program  was  initiat?  i  to 
study  the  failure  dynamics  using  intentionally  induced  failur  This 
report  covers  only  the  surface  temperature  history  of  a  drill?,  hole 
as  a  function  of  time  and  propellant  gas  flow.  An  analytical  -■  del 
was  constructed  frcm  existing  solutions  to  previously  well  b  v.isti- 
gated  phenomena. 

The  utilization  of  this  model  predicts  that  the  surface  of  a  0.  0135 
inch  hole  in  an  aluminum  case  reaches  the  melting  point  ear’.y  in  the 
internal  ballistic  cycle  when  otherwise  unrestricted  flow  oi  propellant 
combustion  gas  and  products  through  the  hole  are  allowed  to  ensue. 

Application  of  this  same  analysis  to  a  brass  case  indicates  that 
the  melting  point  of  the  surface  of  a  hole  in  the  case  is  only  approached 
at  a  point  in  time  late  in  the  internal  ballistic  cycle. 
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vThe  consequence  of  propellant  gas  flow  through  a  split  in  the  cartridge  case 
sidewall  is  usually  more  serious  with  aluminum  cases  than  it  is  with  brass  cases. 
In  order  to  gain  a  better  understanding  of  the  failure  phenomenon  experienced  witl 
aluminum  cases,  a  combination  experimental  and  theoretical  program  was  initiate 
ed  to  study  the  failure  dynamics  using  intentionally  induced  failures.  This  report 
covers  only  the  surface  temperature  history  of  a  drilled  hole  as  function  of  time 
and  propellant  gas  flow.  An  analytical  model  was  constructed  from  existing  solu¬ 
tions  to  previously  well  investigated  phenomena. 

The  utilization  of  this  model  predicts  that  the  surface  of  a  0.  0135  inch  hole  in 
an  aluminum  case  reaches  the  melting  point  early  in  the  internal  ballistic  cycle 
when  otherwise  unrestricted  flow  of  propellant  combustion  gas  and  products 
through  the  hole  are  allowed  to  ensue. 

Application  of  this  same  analysis  to  a  brass  case  indicates  that  the  melting  pai: 
of  the  surface  of  a  hole  in  the  case  is  only  approached  at  a  point  in  time  late  in  th< 
internal  ballistic  cycle. 
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GLOSSARY 


A  Surface  area  of  bore 

Ag  Surface  area  of  radiating  gas  column 
C  Initial  propellant  charge 

Cp  Specific  heat  at  constant  pressure 

Cv  Specific  heat  at  constant  volume 

D  Diameter  of  bore 

E  Constant  used  in  Reynolds  analogy 

K  Stefan-Boltzmann  constant 

L  Length  of  bore 

M  Local  Mach  number 

Mach  number  at  position  i 
MW  Molecular  weight  of  propellant  gases 
Nuz  Nusselt  number 

P  Gas  pressure 

Pr  Prandtl  number 

Q  Radiated  energy 

Qgas-*-solid  Total  thermal  energy  transferred  from  the  gas  to  the  solid 
R  Gas  constant 

Re 2  Reynolds  number  (based  on  lengch) 

RF  Recovery  factor 

T  Gas  temperature 

Taw  Adiabatic  wall  temperature 

Initial  temperature  of  solid 
Tg  Characteristic  gas  temperature 

Ts  Temperature  in  solid 

T#  Reference  temperature 

Change  in  temperature  due  to  radiation 

a  Radius  of  bore 

f  Friction  factor 

hcz  Heat  transfer  coefficient 

k  Thermal  conductivity 
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GLOSSARY  -  Cont'd 


kg  Gas  thermal  conductivity 

i  Length  dimension  defined  by  v  (  At) 

n  Number  of  moles  of  propellant  gas 
q  Heat  flux 

r  Distance  in  radial  direction  from  bore's  centerline 

t  Time 

^t  Time  interval  for  radiation  to  occur 

v  Gas  velocity 

x  Distance  into  solid  (cartesian  variable) 

z  Distance  along  bore's  axis 

5  Boundary  layer  thickness 

Boundary  layer  displacement  thickness 
e(t)  Thickness  of  thermal  layer 

y  Ratio  of  specific  heats 

K  Thermal  diffusi/ity  in  solid 

p  Gas  density 

p  Absolute  viscosity  of  propellent  gas 

Subscripts 

o  Stagnation  state  (Reservoir) 

s  Properties  characteristic  of  the  solid 

oo  Free  stream  condition 


INTRODUCTION 


The  United  States  Army  is  performing  exploratory  development 
to  determine  the  engineering  parameters  required  for  reliable  use  of 
aluminum  cartridge  cases  in  high  pressure,  small  caliber-ammunition 
applications.  This  work  is  based  on  the  need  to  create  lightweight 
weapon  systems  in  near  future  rpolications.  Conservation  of  copper 
resources  is  an  added  benefit  to  be  gained  from  success  in  this  effort. 

Past  experience  has  shown  that  direct  substitution  of  aluminum 
cased  ammunition  in  existing  weapon  systems  po|es  a  potential  prob¬ 
lem  with  respect  to  the  case  failure  mechanism.  In  certain  .instances 
a  case  failure,  commonly  called  "burn -through,  "  results  in  gas-washed 
and  eroded  case  surface  and  represents  a  possible  damage  mechanism 
to  both  weapon  operator  and  weapon.  Evidence  to  date  suggests  that 
this  phenomenon  results  from  a  mechanical  failure  of  the  cartridge 
case  during  the  internal  ballistic  cycle  or  firing  of  the  weapon. 

When  a  brass  case  fails  during  firing,  the  rifleman  often  is  una¬ 
ware  that  a  failure  has  occurred.  The  usual  consequences  of  such  a 
failure  are  slight  gas  leakage  at  the  weapon  breech,  discoloration  of  the 
case  at  the  point  of  the  defect  and,  possibly,  a  slight  change  in  cyclic 
rate  of  the  weapon  when  such  firing  is  in  the  automatic  mode.  On  the 
other  hand,  the  firing  of  ammunition  having  a  defective  aluminum  car¬ 
tridge  case  can  be  far  more  spectacular.  Here,  there  can  be  an 
efflux  of  very  luminous  gases  at  the  weapon  breech  which  results  in 
weapon  and  case  damage;  the  rifleman  is  also  exposed  to  hazard  under 
these  circumstances. 

The  positive  approach  taken  to  serve  as  the  basis  for  successful 
application  of  aluminum  cases  consists  in  understanding  and  solving  the 
problem.  Prior  to  initiation  of  the  "understanding  the  problem"  phase, 
various  suppositions  were  expressed  as  to  the  mechanism  of  this  phe¬ 
nomenon.  One  technique  was  found  to  be  reliable  in  producing  the 
characteristic  flash  and  metal  erosion  which  has  sometimes  been  ob¬ 
served  during  firing  of  defective  aluminum  cased  ammunition. 

^Miller,  S. ,  "Design,  Development  and  Fabrication  of  100,000  Car¬ 
tridges,  Ball,  Caliber  .50,  M33  Type  Assembled  with  Case,  Cartridge 
Aluminum  .Caliber  .50,  FAT  39,  "  Technical  Report  R-1265,  June  1956, 
Frankford  Arsenal,  Philadelphia,  Pa. 
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Summerfield,  M.  "Letter  to  Reed  E.  DonnardjSubject;  'Cartridge  Burn 
Through  Problem',  "  Sept  1969,  Princeton  University,  Guggenheim  Aero 
space  Propulsion  Lab.  ,  Princeton,  New  Jersey. 
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This  technique  consisted  in  preplacing  a  gas  path  in  the  case  by  drill¬ 
ing  a  hole  through  the  extractor  groove  region  to  the  interior  of  the 
case  before  the  case  was  assembled  into  a  cartridge.  Thus,  this  in¬ 
tentional  case  defect  served  as  an  excellent  means  for  observing  this 
phenomenon  called  "burn-through"  in  the  actual  weapon /ammunition 
environment.  For,  after  all,  this  represented  the  "real  world"  in 
which  the  problem  existed  and  in  which  it  had  to  be  solved.  Aluminum 
cased  cartridges  fired  in  this  manner  exhibited  the  characteristics 
normally  identified  as  a  'burn-through.  "  Other  techniques  were  used 
to  study  this  failure  mechanism  but  they  will  be  discussed  in  another 
report. 

An  observation  that  propellant  gas  flow  through  a  path  in  the  case 
is  the  underlying  process  whereby  "burn -through"  is  manifested  was 
demonstrated  in  an  experiment  wherein  propellant  gas  flow  through  an 
induced  path  in  the  case  was  reduced  to  almost  nil  by  the  use  of  a  spe¬ 
cial  weapon  bolt/ chamber  design.  There  was  no  damage  to  the  case 
or  weapon  test  fixture  even  though  a  gas  path  existed  in  the  case  wall 
at  a  location  where  typical  case  and  weapon  damage  would  have  other¬ 
wise  been  obtained.  Stopping  or,  at  least,  drastically  reducing  pro¬ 
pellant  gas  flow  through  the  path  in  the  case  eliminated  the  damage,  ^ 
Since  the  physical  properties  of  aluminum  and  brass  differ  markedly 
(density,  specific  heat  and  melting  point),  the  results  of  propellant  gas 
flow  through  a  path  in  cases  made  from  these  materials  also  differ  sig¬ 
nificantly.  Table  I  comparing  these  properties  was  abstracted  from 
reference  4.  In  addition  to  these  property  differences,  chemical  activ¬ 
ity  of  aluminum  alloy  is  much  greater  than  that  of  brass. 


The  study  to  identify  the  parameters  associated  with  the  failure 
mechanism  of  the  aluminum  case  and  to  ascertain  cheir  effect  on  the 
severity  of  resultant  system  damage  is,  by  necessity,  broad  in  nature. 
Such  aspects  as  the  relationship  of  propellant  gas  flow,  energy  trans¬ 
fer  from  the  propellant  gas  to  exposed  case  surface,  heating  response 
of  the  exposed  case  surface  and  chemical  reactivity  of  the  aluminum  case 
material  are  under  investigation. 

The  subject  matter  of  this  report  is  concerned  with  one  part  or 
building  block  in  the  model  that  will  evolve  to  describe  this  phenomenon. 
Material  contained  herein  describes  the  local  or  surface  temperature 
profile  along  the  gas  path  (here  a  drilled  hole)  as  a  function  of  time  and 
the  energy  transfer  available  from  the  combusting  propellant. 


^Unpublished  data,  "The  Aluminum  Cartridge  Case  'Burn -Through' 
Problem- Characteristics,  Isolation,  and  Means  of  ^elimination,  " 
Frankford  Arsenal, 

^Lyman,  T. ,  ed.  ,  Metals  Handbook,  Propeic'fs  and  Selection  of  Metals, 
8th  ed. ,  Vol.  1,  American  Society  for  Metals,  Novelty,  Ohio,~T96T, 
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TABLE  I. 

Comparison  of  70-30  Brass  and  7075  Aluminum 


Property 

Density  (lbm/ft^) 

70-30  Brass 

532.52 

7075  Aluminum 

174.80 

Thermal  conductivity 
(BTU  -  ft/ft2  hr  °F) 

70.00 

70. 15 

Specific  heat  (BTU/lbm  °F) 

0.09 

0.23 

Thermal  diffusivity  (ft^/hr) 

1.46 

1.74 

Solidus  temperature  (*F) 

1680 

890 

Liquidus  temperature  (8F) 

1750 

1180 

STATEMENT  OF  THE  PROBLEM 


The  mathematics  developed  herein  describe  the  heat  flux  resulting 
from  the  fast  moving  propellant  gases  through  the  induced  orifice  and 
allow  determination  of  t’  e  temperature  profiles  in  selected  regions 
of  the  cartridge  case,  emphasis  is  placed  on  the  gas /solid  interface, 
the  interior  surface  the  induced  orifice.  Temperatures  at  the  in¬ 
terior  surface  of  the  hole  are  inferred  to  as  "local"  temperatures. 

To  test  the  accuracy  of  the  mathematica'  model,  the  appropriate 
physical  constants  and  values  obtained  from  the  experiment  were 
substituted  into  the  equations;  temperature  profiles  were  also 
determined. 

Throe  th  orios  have  been  postulated  to  account  for  the  failure  of 
the  aluminium  eaitridge  case.  All  of  these  theories  are,  based  on  the 
heating  of  the  case  as  s^+ed  below: 

1.  It  may  be  possible  that  the  failure  is  the  result  of  pure  melti  lg, 

2.  It  may  be  possible  that  the  failure  is  the  result  of  decreased 
yield  strength  as  a  result  of  increased  temperature.  Then,  solid  and 
possibly  vaporized  aluminum  will  be  part  of  the  gas  stream, 

3.  It  may  be  possible  that  chemical  reactions,  either  occurring 
at  the  surface  of  the  cartridge  case  or  in  close  proximity  to  the  case, 
are  initiated.  This  last  possibility  affects  the  heat  flux  to  the  cartridge 
case. 

To  determine  the  validity  oi  the  model  developed  in  this  document, 
uhe  cartridge  case  will  be  assumed  to  remain  intact  (until  melting)  and 
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that  there  are  no  chemical  reactions.  In  the  event  further  experimen¬ 
tation  indicates  that  the  cartridge  case  undergoes  structural  degrada¬ 
tions  -  at  a  temperature  lower  than  the  melting  point  -  the  analysi 9 
contained  herein  will  also  enable  identifier don  of  that  critical  temper¬ 
ature.  Further  assumptions  are: 

1.  Propellant  gas  behaves  as  an  ideal  gas, 

2.  Classical  heat  conduction  equation  describes  the  dissipation  of 
energy  in  the  solid  until  melting, 

3.  The  physical  properties  of  the  cartridge  case  material  are  con¬ 
stant  during  the  experiment, 

4.  It  is  possible  to  identify  an  average  operating  pressure, 

5.  The  gas  temperature  is  the  adiabatic  flame  temperature, 

To  reiterate,  it  is  only  the  object  of  this  study  to  concentrate  on 
the  temperature  profile  aspects  of  the  problem. 


EXPERIMENTAL 

The  failure  dynamics  were  studied  by  inducing  failures.  Figure  1 
shows  a  typical  5.  56  mm  cartridge  case  and  the  location  where  small 
fissures  were  introduced  in  the  head  regions  of  both  aluminum  alloy  and 
brass  cartridge  cases.  Both  types  of  cartridge  case,  with  the  small  hole 
present,  were  loaded  with  a  charge  of  27.  0  grains  of  WC  846  propellant 
and  fired  in  a  test  weapon.  Although  the  size  and  location  of  this  hole 
might  have  been  varied,  this  study  is  addressed  to  the  hole  geometry 
described  in  Figure  1.  The  thermochemical  and  thermodynamic  pro¬ 
perties  of  WC  846  propellant  are  tabulated  in  Appendix  A. 

The  results  of  this  testing  are  shown  in  Figures  2  and  3.  As  seen 
in  Figure  2,  the  two  brass  cases  seem  unaffected  after  firing.  However, 
the  two  aluminum  cartridge  cases  have  noticeable  damage  in  the  head 
region.  The  unfired  aluminum  alloy  cartridge  case  can  be  used  to  com¬ 
pare  the  damage  after  firing  with  its  initial,  drilled  condition.  Figure 
3  shows  sectioned  aluminum  alloy  and  brass  cartridge  cases  before  and 
after  firir  .  As  indicated  in  the  photograph,  the  small  fissure  in  the 
brass  case  undergoes  only  a  slight  alteration.  On  the  other  hand,  the 
increase  in  the  hole's  diameter  of  the  aluminum  alloy  cartridge  case  is 
rather  obvious. 

After  initiation  of  the  primer  by  the  weapon's  firing  pin,  the  pro¬ 
pellant  grains  burn,  releasing  gas  and  increasing  the  pressure.  This 
pressure  is  of  sufficient  magnitude  to  force  the  propellant  gases  through 
any  unsealed  opening  in  the  cartridge  case.  When  this  opening  is  in  the 
head  region  of  the  cartridge  case,  the  gas  flow  is  as  indicated  in  Figure  4. 
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The  generation  of  gas  has  the  additional  effect  of  expanding  the 
sidewalls  of  the  case  until  they  physically  touch  the  weapon's  chamber. 
Although  modern  day  weapons  are  designed  so  that  the  chamber  offers 
the  maximum  possible  support  to  the  cartridge  case,  there  is  a  small 
region  which  is  totally  unsupported.  The  existence  of  this  region  is  a 
requirement  imposed  on  the  cartridge  case  by  the  weapon  designer  to 
allow  insertion  of  the  case  into  the  chamber  and  extraction  after  firing. 
It  has  been  observed  that  the  erosion  of  case  material  and  damage  to 
the  weapon  is  most  severe  in  this  region  where  the  cartridge  case  is 
unsupported  by  the  gun  chamber.  This  elimination  of  gas  flow  and  sub¬ 
sequent  observation  of  minimal  damage  is  borne  out  in  the  study  where¬ 
in  the  weapon  was  modified  to  give  total  support  to  the  cartridge  case. 
That  is,  the  automatic  firing  mode  of  the  weapon  was  sacrified  by  using 
the  weapon's  bolt  to  surround  and  support  entirely  the  region  of  the  case 
reserved  for  chambering  and  extraction. 


To  obtain  pressure  versus  time  data,  a  thick-walled  test  barrel 
was  drilled  and  threaded  to  accept  a  607A  Kistler  gauge.  This  gauge 
was  positioned  midway  along  the  longitudinal  axis  of  the  cartridge  case 
to  measure  the  gas  pressure.  Figure  5  shows  the  pressure-time 
history  of  the  propellant  gases  inside  a  5.  56  mm  cartridge  case  which 
contained  the  charge  of  27.0  grains  of  WC  846  propellant.  As  can 
readily  be  seen  the  pressure  increases  from  0  to  50,  000  psi  in 
approximately  0.6  milliseconds  and  the  entire  event  is  known  to 
last  for  approximately  2.0  milliseconds. 

In  addition  to  the  pressure-time  curve,  it  was  desired  to  monitor 
the  duration  of  the  gas  flow  from  the  orifice.  This  data  was  obtained 
by  machining  a  test  weapon  to  allow  exposure  of  the  head  region  of  the 
cartridge  case.  Previous  visual  observation  of  the  failure  process 
showed  the  gases  which  exit  the  induced  orifice  in  the  cartridge  case 
were  luminous.  Hence,  positioning  a  photoelectric  cell  in  the  same 
horizontal  plane  and  at  right  angles  to  the  induced  orifice  provided 
the  desired  information.  That  is,  by  recording  the  output  of  the  photo¬ 
electric  cell  on  the  screen  of  an  oscilloscope  it  is  possible  to  determine 
when  the  gases  first  exit  the  orifice  and  for  how  long  the  flow  continues. 
A  dual  beam  oscilloscope  was  used  to  record  the  pressure-time  and 
the  photoelectric  cell  output-time  histories  simultaneously.  Analysis 
of  this  data  indicates  that  gas  flow  begins  with  the  first  significant 
pressure  rise  and  continues  until  the  pressure  curve  returns  to  zero 
pressure.  Hence,  there  is  gas  flowing  from  the  orifice  for  a  period 
of  2.0  milliseconds. 


O 

JUnpublished  data,  "The  Aluminum  Cartridge  Case  'Burn-Through' 
Problem  -  Characteristics,  Isolation,  and  Means  of  Elimination,  " 
Frankford  Arsenal,  Philadelphia,  Penna. 
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THEORETICAL  ANALYSIS 


The  Flow 

The  approach  used  was  to  develop  a  mathematical  model  and  to 
calculate  the  heat  transfer  rate  and  temperature  increases  at  the 
orifice  walls  to  determine  if  melLng  temperatures  can  be  reached. 
Hence,  the  analysis  which  follows  is  directed  to  the  cartridge  case 
as  shown  in  Figure  1.  As  is  seen  in  this  figure,  the  hole  is  drilled 
in  the  most  massive  region  of  the  case.  Also,  the  L/D  ratio  of  the 
hole  is  approximately  5.  Therefore,  a  reasonable  approximation 
to  the  region  of  the  cas'j  where  the  hole  is  drilled  is  that  of  a  small 
hole  in  a  medium  of  infinite  extent  in  the  radial  direction.  This 
description  is  best  seen  by  referring  to  Figure  6.  Distance  along 
the  axis  of  the  bore  is  specified  by  the  variable  z,  and  distance 
from  the  center  line  of  the  bore  in  the  radial  direction  is  given  by 
the  variable  r. 

A  knowledge  of  the  propellant  gas  pressure  and  temperature, 
coupled  with  the  appropriate  flow  theory  in  the  bore,  leads  to  a 
description  of  the  heat  flux  incident  on  the  bore1  s  surface.  The 
classical  heat  diffusion  equation,  written  in  the  appropriate  coor¬ 
dinate  system,  can  be  solved  subject  to  boundary  conditions  de¬ 
termined  from  the  flow  analysis. 

It  is  first  necessary  to  describe  the  ilow  of  propellant  gases 
from  the  interior  of  the  cartridge  case  through  the  induced  fissure. 

To  do  so,  it  is  important  to  select  the  most  representative  gas  dynamic 
model  based  upon  bore  size.  With  a  L/D  ratio  of  approximately  5, 
the  relatively  small  opening  of  0.  0135  inch  in  diameter,  and  the 
rapid  rate  of  pressurization,  the  flow  may  be  assumed  to  be  choked. 

The  additional  assumption  is  made  that  the  flow  is  developed  to  the 
point  where  it  consists  of  a  well  defined  boundary  layer  and  a  central, 
core  region.  Figure  7  is  a  drawing  depicting  this  flow  process.  The 
temperature,  pressure,  density,  and  velocity  as  a  function  of  position 
in  the  core  region  can  be  determined  with  the  appropriate  flow  model. 
Determination  of  the  nature  of  the  boundary  layer,  laminar  or  turbulent, 
can  be  accomplished  by  an  examination  of  the  magnitude  of  the  local 
Reynolds  number. 

Lee  and  Sears^  suggest  an  adiabatic  treatment  of  the  flow  be 
considered  providing  the  bore  is  not  too  long  (less  than  ten  bore 
diameters).  Admittedly,  the  flow  process  should  a  ccount  for  the 
substantial  transfer  of  energy  to  the  bore's  sidewalls  since  melting 
has  been  hypothesized  to  occur.  The  adiabatic  assumption  of  the 


5Lee,  J.  F.,  and  Sears,  F.  W.,  Thermodynamics,  Addison-Wesley, 
Cambridge,  Mass,,  1955. 
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Schematic 


flow  does  not  contradict  the  hypothesized  fact  that  there  is  sub¬ 
stantial  energy  transfer.  The  analysis  of  the  core  region  merely 
permits  a  determination  of  the  gas  conditions  -  i.  e. ,  temperature, 
pressure,  density,  and  velocity  -  external  to  the  boundary  layer. 

It  must  be  remembered  ihat  it  is  across  the  boundary  layer  that  the 
energy  is  transferred  to  the  bore's  sidewalls.  The  gas  dynamics  in 
the  core  will  be  treated  as  adiabatic  flow  with  friction  in  a  duct  of 
constant  area. 

£ 

Shapiro  developed  a  series  of  working  formulas  to  describe 
adiabatic  flow  with  friction  in  a  duct  of  constant  area.  Momentum, 
energy,  and  mass  equations  are  written  for  the  flow  of  a  perfect  gas 
through  an  elemental  control  volume.  Since  the  reservoir  (chamber) 
conditions  are  available  and  since  choking  is  assumed  to  occur  at  the 
exit,  these  formulas  may  be  used  to  predict  the  temperature,  pressure, 
velocity,  and  density  at  any  position  along  the  bore's  axis.  A  second 
alternative  is  available.  Keenan  and  Kaye?  provide  a  series  of  tables 
for  one -dimensional  adiabatic  flow  with  friction  through  a  constant 
area  duct.  Although  these  tables  are  available  for  gases  of  differing 
ratio  of  specific  heats,  they  are  all  based  on  an  Ideal  Gas  Equation 
of  State.  An  Ideal  Gas  Equation  of  State  is  used  throughout  this 
analysis  to  relate  the  intrinsic  parameters  of  the  system. 

The  Mach  number  at  any  position  along  the  bore's  .  xis  is  now  to 
be  determined.  Lee  and  Sears5  give  a  relationship  for  .he  bore  length 
L,  reauired  for  the  flow  to  pass  from  a  Mach  number  Mj,  to  a  Mach 
number  M2  as 

L  =  (Lmax)Mj  -  (Lnnax)M2  . 

Multiplication  bv  4f/D  on  both  sides  of  the  equation  changes  the 
form  to 


_  4f 

D  D  '  maxi 

where  f  is  the  friction  factor 

D  is  the  diameter  of  the  bore, 
obtained  from  the  Reynolds  analogy 

‘Vz  _  F  f 
pvCp  -  E  2 

o - : —  1  ■ 1  ■ 1  ~ 

Lee,  J,  F.  ,  and  Sears,  F,  W.,  Thermodynamics,  Addison-Wesley, 
Cambridge,  Mass.,  1955. 

^Shapiro,  A.  H. ,  The  Dynamics  and  Thermodynamics  of  Compressible 
Fluid  Flow,  Vol.  1,  Ronald  Press,  New  York,  1953. 

?K  eenan,  J.  H.  ,  and  Kaye,  J.  ,  Gas  Tables,  Wiley,  New  York,  1956 
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where  hcz  is  the  film  coefficient  of  heat  transfer  or  simply 
the  heat  transfer  coefficient 


p  is  the  gas  density 
v  is  the  gas  velocity 

Cp  is  the  specific  heat  at  constant  pressure. 
With  E  =  1  the  above  equation  can  be  rearranged  to  read 


f 


^  ^cz 

PvCp 


(2) 


However,  the  friction  factor  is  unknown  at  this  time  in  the  analysis. 
Shapiro®  states  friction  factors  in  the  range  of  0.001  to  0.004  are 
realistic  for  the  type  of  flow  with  which  this  analysis  is  concerned. 
Keenan  and  Kaye's  '  gas  tables  are  used  to  determine  the  Mach 
number  as  a  function  of  position  along  the  bore  for  friction  factors 
of  0.001,  0.002,  and  0.004.  For  the  time  being,  these  three  values 
of  the  friction  factor  will  be  used  since  hcz,  p,  and  v  are  all  unknown. 
Later  in  the  analysis  this  requirement  will  be  relaxed  and  only  one 
value  will  be  used.  The  ratio  of  specific  heats  for  propellant  gas 
given  in  Appendix  A  is  1. 24.  Since  Keenan  and  Kaye7  do  not  provide 
a  separate  table  for  y=  1.  24,  it  is  necessary  to  interpolate  between 
y =  1. 20  and  y  -  1. 30. 

With  M2  =  1,  the  ratio  (4f/D)  (Lmax)M2  becomes  0  and  equation 
(1)  reduces  to 


if  L 
D 


4f 

D 


(LmaxJjvty 


Values  of  4fL/D  for  f  =  0.  001,  0.  002,  and  0.  004,  L  =  0.  062  inch  and 
D  =  0.0135,  are  given  in  Table  II.  The  corresponding  Mach  numbers 
for  the  particular  friction  factor  are  obtained  from  the  body  of  the  gas 
tables.  Since  the  length  used  to  calculate  4fL/D  was  0.062  inch  (the 
bore's  length),  these  Mach  numbers  are  the  inlet  Mach  numbers. 


Equation  (1)  is  used  again  to  calculate  the  distance  required  for 
the  Mach  number  to  change  from  the  inlet  value  to  some  intermediate 
value.  A  typical  analysis  for  a  friction  factor  of  f  -  0.  002  is  given 
in  Appendix  B.  Figure  8  shows  the  dependence  of  local  Mach  number 
on  the  position  along  the  longitudinal  axis  of  the  bore  for  the  three 
friction  factors. 


®  3h*  oiro,  A.  H. ,  The  Dynamics  and  Thermodynamics  of  Compressible 
Fluid  Flow,  Vol.  1,  Ronald  Press,  New  York,  1953. 

Keenan,  J.  H.  and  Kaye,  J.  ,  Gas  Tables,  Wiley,  New  York,  1956. 
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POSITION  ALONG  BORE'S  AXIS  (INCHES) 

Figure  8,  Mach  Number  as  a  Function  of  Position  Along  the  Bore 
(F  fiction  Factor  as  Indicated) 
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TABLE  n. 

Inlet  Mach  Numbers 


Friction  Factor  4fL/D  Mach  Number 

0.001  0.0183  0.900 

0.002  '0.0367  0.865 

0.004  0.0734  0.815 


Once  the  local  Mach  number  is  available,  it  is  novr  possible  to 
evaluate  both  the  free  stream  and  stagnation  values  for  the  tempera¬ 
ture,  pressure,-  and  density  at  any  point  along  the  bore.  The  chamber 
conditions,  as  given  in  Appendix  A,  to  be  used  throughout  this  analysis 
are: 


T0  =  5040  °R 
P0  =  2.  5  x  10^  psi 
R  =  64.  372  ft-ibf/lbm°R 
y  =1.24 

The  subscript  o  refers  to  a  stagnate  n  state.  With  an  Ideal  Gas 
Equation  of  State,  the  gas  density  in  u.  3  chamber  (  stagnation  state  ) 
is  given  by 

Po  lbm 

Po  =  It;  =  11.09  -jjf*  •  P> 

Isentropic  conditions  are  thought  to  be  valid  at  the  inlet  of  the  bore. 
The  isentropic  conditions  merely  account  for  the  flow  of  propellant  gas 
from  the  chamber  to  the  beginning  of  the  bore  -  the  inlet.  For  the 
isentropic  flow,  Shapiro ^  provides  the  following  governing  equations; 


(4) 


Poo  = 


1  +  M2 


y/v-1 


(5) 


^Shapiro,  A.  H. ,  The  Dynamics  and  Thermc  dynamics  of  Compressible 
Fluid  Flow,  Vol.  1,  Ronald  Press,  New  York,  1953. 
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Pco 


(6) 


_ Po_ 

(1+-2£lM2) 

Free  stream  conditions  are  identified  v/ith  the  subscript  “  .  The  local 
value  for  the  gas  velocity  may  also  be  determined  from  the  formula 


v  =  M  y/yRT ~  .  (7) 

Using  equations  (4),  (5),  (6),  and  (7),  with  T0,  PQ,  and  pQ  as 
previously  specified  and  the  inlet  Mach  number  obtained  from  Figure  8, 
it  is  possible  to  specify  completely  all  the  conditions  at  the  inlet  of  the 
bore.  This  information  is  presented  in  Table  III. 


TABLE  III. 

Gas  Properties  at  the  Inlet  Assuming  Isentzopic  Process 


Flic t ion 

Mach 

Stagnation 

Free  Stream 

Gas 

Factor 

Number 

Tempt  tture 

Pressure 

Density 

Temperature 

Pressure  Density 

Velocity 

(-) 

(-) 

(°bj 

(psi) 

(’bn/ ft') 

(°M 

(P»i)  (lbm/ft3) 

(ft/sec) 

0.001 

0.90 

5.040  x  103 

2.5  x  10* 

11.09 

4.593  x 

103 

1.524  x  104  7.42 

3.095  x  103 

0.002 

0.865 

5.040  x  103 

2.5  x  10* 

11.09 

4.624  x 

103 

1.602  x  104  7.75 

2.981  x  103 

0.004 

0.815 

5.040  x  103 

2.5  x  104 

11.09 

4.667  x 

103 

1.677  x  104  8.04 

2.822  x  103 

Attention  is  now  turned  to  the  value  of  the  oarameters,  identified 
in  Table  III,  at  various  locations  in  the  bore.  Characteristic  of  the 
adiabatic,  flow  is  the  fact  that  the  stagnation  temperature  remains 
constant  throughout  the  bore.  Equations  (4),  (5),  (6),  and  (7)  and  the 
gas  tables  are  used  to  determine  the  gas  properties  at  the  intermediate 
locations.  These  calculations  are  performed  for  a  friction  factor  of 
0.  002  at  positions  where  the  Mach  numbers  are  0.  90,  0.  95  and  1.  00 
ard  are  given  in  Appendix  C.  This  data  is  shown  in  Figure  9. 

At  this  point  in  the  analysis,  an  approximation  is  used.  Basically, 
the  problem  which  is  being  studied  is  that  of  flow  and  heat  transfer  in 
a  small  bore.  Kreith®  states  that  for  very  short  tubes  or  rectangular 
ducts  with  initially  uniform  velocity  and  temperature  distribution,  the 
flow  conditions  along  the  wall  approximate  those  along  a  flat  plate. 


®Kreith,  F.,  Principles  of  Heat  Transfer,  International  Textbook 
Company,  Scranton,  Penna. ,  1963. 
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Hence,  the  original  two  dimensional  cylindrical  geometry  (r,  z)  can 
be  replaced  with  a  two  dimensional  cartesian  system  (x,  z).  Further 
justification  for  the  flat  plate  treatment  of  the  problem  may  be  obtained 
by  showing  that  the  boundary  layer  displacement  thickness  is  small 
compared  to  the  bore's  radius.  This  fact  implies  that  the  boundary 
layer  is  essentially  localized  near  the  surface  of  the  bore. 

Since  the  gas  properties  have  been  determined  (Figure  9).  it  is 
possible  to  evaluate  the  Prandtl  and  Reynolds  numbers.  The  values 
for  the  specific  heat  at  constant  pressure,  absolute  viscosity,  and 
thermal  conductivity  are  available  in  Appendix  A.  The  Prandtl 
number,  defined  by  Pr  =  Cpp/kg  has  a  value  of  2.  905.  To  de¬ 
termine  the  nature  of  the  boundary  layer,  turbulent  or  laminar, 
the  local  Reynolds  number  is  needed.  The  Reynolds  number,  based 
on  length,  is  given  by  Rez  =  vpz/p.  The  values  of  the  velocity,  v, 
and  the  gas  density,  p ,  for  a  particular  value  of  z  are  obtained  from 
Figure  9.  For  example,  at  z  =  0.031  inch  (the  bore  s  midpoint)  the 
local  Reynolds  number  is  1.26  x  10&;  likewise  at  z  =  0.062  inch  (the 
bore's  exit)  the  local  Reynolds  number  is  2.40  x  10°.  Kreith®  reports 
that  the  flow  over  a  flat  plate  is  turbulent  where  the  local  Reynolds 
number  is  approximately  3  x  10^,  Comparing  3  x  10^  with  the  cal¬ 
culated  values  of  1.26  x  10^  (z  =  0.031  inch)  and  2.40  x  10°  (z  =  0.062 
inch),  it  is  possible  to  conclude  that  the  flow  is  indeed  turbulent. 

It  is  now  possible  to  calculate  the  boundary  layer  thickness  at 
any  point  along  the  axis  of  the  bore.  Rohsenow  and  Choi9  report 
that  the  boundary  layer  thickness  for  turbulent  flow  is  given  by 


6 


0.  37  z 


(Rez) 


1/5 


(8) 


For  the  flat  plate  assumption  to  be  valid  it  is  necessary  to  show  that 


a 


(9) 


where  a  is  the  bore's  radius  and  6*1?  the  boundary  layer  displace¬ 
ment  thickness.  At  z  =  0.031  inch,  the  local  Reynolds  number,  Rez, 
is  1.26  x  106.  Olson8 * 10  states  that  the  boundary  layer  displacement 
thickness  for  turbulent  flow  is  given  by  6/8;  hence,  equations  (8)  and 
(9)  may  be  combined  to  read 


8Kreith,  F.,  Principles  of  Heat  Transfer,  International  Textbook 
Company,  Sc ranton,  Penna. ,  1 9^3. 

9Rohsenow,  W.  M. ,  and  Choi,  H. ,  Heat,  Mass,  and  Momentum 
Transfer,  Prentice-Hall,  Englewood  Cliffs,  N.  J. ,  19&1. 

10Olson,  R.  M. ,  Essentials  of  Engineering  Fluid  Mechanics,  Inter 
national  Textbook  Company,  Scranton,  Penna.,  1962. 
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6*  = 


(10) 


0.  37  z 

1  /5  • 

8  a  (Reg.)1  ^ 

Evaluating  equation  (10)  at  z  =  0.031  inch  and  a  =  0.00675  inch,  the 
following  result  is  obtained:  6*/a  ?  0.0128.  Since  6^/a  has  been  de¬ 
termined  to  be  0.  0128,  which  is  very  much  less  than  one,  the  flat 
plate  treatment  of  the  bore's  surface  is  justified. 

A  great  deal  of  attention  has  been  given  to  the  flow  character¬ 
istics  of  the  propellant  gases  in  the  small  bore.  This  is  attributable 
to  the  fact  that  the  flow  has  been  identified  as  being  the  means  whereby 
the  energy  is  transferred  to  the  bore's  sidewalls.  This  energy  trans¬ 
port  can  take  place  by  basically  two  means,  forced  convecticn  and 
radiation.  It  is  therefore  important  to  determine  the  Nusselt  number 
(the  dimensionless  quantity  relating  heat  transfer  and  fluid  flow)  in 
terms  of  the  Prandtl  and  Reynolds  numbers.  With  turbulent  flow 
over  a  flat  plate,  OlsonlO  states  that  a  reasonable  value  for  the 
Nusselt  number  may  be  obtained  from 

Nuz  =  0.  0288  (Pr)l/3  (Rez)4/5  (11) 

where  Nuz  implies  that  the  Nusselt  number  is  a  function  of  position 
along  the  bore's  longitudinal  axis. 

With  both  the  Prandtl  and  Reynolds  numbers  known  at  this  time 
and  with  the  classical  definition  of  the  Nusselt  number  as 


Nuz  = 


where 


is  the  heat  transfer  coefficient 

is  the  thermal  conductivity  of  the  propellant  gas 


(12) 


it  is  possible  to  check  on  the  internal  consistency  of  the  friction 
factor,  f  .  Substitution  of  equation  (11)  into  equation  (12)  and  re¬ 
arrangement  permits  identification  of  hcz  in  terms  of  known  quantities 


hcz  =  (0.  0288  )(Pr)1/3  (Rez)4/5 

Equation  (2)  now  reads 

f  -  (2)  (0-  0288)  (Pr)1/3(Rez)4/5  kg 

PC  pvz 


(13) 


^Olson,  R.  M. ,  Essentials  of  Engineering  Fluid  Mechanics,  Inter¬ 
national  Textbook  Company,  Scranton,  Penna.  ,  1 962. 
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Evaluating  equation  (13)  at  the  bore's  midpoint,  z  =  0.031  inch,  and 
obtaining  the  required  values  from  Figure  9  and  Appendix  A  a  realistic 
friction  factor  is  found  to  be  f  =  0.  002. 

A  review  of  the  analysis  thus  far  is  given.  This  section  has 
addressed  to  the  adiabatic  flow  of  a  perfect  gas  through  a  small  bore 
choked  at  the  exit  plane.  A  characteristic  design  length  equation  was 
used  to  calculate  the  Mach  number  at  any  position  along  the  bore's 
axis.  To  determine  values  for  the  gas  properties  at  the  inlet,  an 
isentropic  flow  process  was  used.  Knowledge  of  the  local  Mach 
number  and  the  inlet  conditions,  assuming  adiabatic  flow  through  a 
duct  of  constant  area  with  friction,  allowed  determination  of  the  gas 
properties  at  any  point  along  the  bore's  axis.  Next,  examination  of 
the  local  Reynolds  number  identified  the  flow  as  being  turbulent. 
Finally,  a  flat  plate  treatment  of  the  bore's  surface  was  justified 
and  the  appropriate  Nusselt  number  constructed. 


The  Heat  Transfer 


As  identified  in  the  previous  section,  the  energy  transport  to  the 
bore's  sidewalls  is  accomplished  by  forced  convection  and  radiation. 

The  heat  which  is  transferred  within  the  solid  is  governed  by  conduction 
and  will  be  treated  later  in  the  analysis.  Since  the  propellant  gases 
are  at  a  relatively  high  temperature,  2800  0  K,  the  effect  of  energy 
transfer  due  to  radiation  cannot  be  disregarded  without  some  justi¬ 
fication.  Therefore,  let  us  calculate  the  radiation  heat  loss  with  the 
assumption  that  the  gas  radiates  like  a  black  body  under  the  worst 
possible  conditions.  Namely,  the  temperature  of  the  gas  as  it  flows 
throughout  the  bore  is  the  adiabatic  flame  temperature,  and  all  the 
propellant  gas  generated  in  the  reaction  flows  through  the  induced 
hole.  (A  closed  bomb  experiment  wherein  the  gas  is  allowed  to  flow 
through  a  hole  0.  0135  inch  in  diameter  for  2.  0  milliseconds. ) 

From  Appendix  A,  the  propellant  properties  required  for  this 
calculation  are: 

Adiabatic  flame  temperature,  TQ  =  5040  °R 

Moles  of  gas  produced,  n  =  0.  0423  — an> — m°  y-r — - 
°  c  gram  of  propellant 

Specific  heat  at  constant  volume,  Cv  =  0.  344  BTU/lbm°R 

Initial  propellant  charge,  C  =  27  grains 

Molecular  weight  of  propellant  gas,  MW  =  23.  62. 

In  general,  with  the  black  body  assumption,  the  radiated  energy,  Q,  is 
Q  =  KTq4  (At)  As  (14) 


22 


where  K 


is  Stefan- Boltzmann  constant  [0. 173  x  10"®  — - — - 

\  ft2-hr(°R)4; 

At  is  time  interval  for  radiation  to  occur 

As  is  surface  area  of  radiating  gas  column. 

The  surface  area  of  the  radiating  gas  column  may  be  calculated  as 
Ag  =  2na  &  with  l  a  length  dimension,  is  defined  as  l  =  v,,,  (At)  and 
where  is  the  greatest  value  for  the  gas  velocity.  The  value  for 
V®  is  obtained  from  Figure  9  and  is  va  =  3400  ft/sec.  Therefore, 
the  surface  area  of  the  radiating  gas  column  is  As  =  2.4  x  10~2  ft^. 
Substituting  this  value  into  equation  (14),  the  radiated  energy  is 
Q  =  202.95  x  10-4  BTU. 


The  radiated  energy  transferred  from  the  propellant  gas  may  be 
set  equal  to  the  internal  energy  change  of  the  gas,  or  Q  =  nCy  AT 
where  AT  is  the  change  in  temperature  as  a  result  of  radiation. 

Since  the  combustion  of  27.  0  grains  of  WC  846  ball  propellant  produces 
0.0738  moles  of  gas,  the  weight  of  this  amount  of  gas  is  given  as  (moles 
of  gas)  x  (molecular  weight)  x  (0.  0022  pounds /gram)  or  0.  0038  pounds. 
The  change  in  temperature  can  now  be  calculated  as 


AT 


Q 

n  Cv 


AT  =  15.49  °F 


Hence,  the  temperature  change  in  the  gas  as  a  result  of  radiation 
under  the  worst  possible  set  of  conditions  is  only  15.49°F.  It  is  there¬ 
fore  concluded  that  the  effect  of  radiation  is  insignificant. 

The  heating-up  of  the  bore  surface  is  accomplished  for  the  most 
part  by  forced  convection  across  a  turbulent  boundary  layer.  Assuming 
that  the  heating-up  varies  only  along  the  axis  of  the  bore,  it  is  possible 
to  write  in  general  terms  that 


Q 


gas¬ 


's  olid 


•A  hcz  (Ts  - 


T  ) 


(15) 


where  Qga.s-*-solid 
A 

^cz 

Ts 


is  the  total  thermal  energy  transferred  from  the 
gas  to  the  solid 

is  the  area  of  the  bore's  surface 
is  the  local  heat  transfer  coefficient 
is  the  temperature  in  the  solid 
is  the  characteristic  gas  temperature  . 
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Olson  points  out  that  the  parameter  hcz  is  extremely  sensitive  and 
depends  on  the  gas  property  values,  such  as  the  specific  heat,  thermal 
conductivity,  viscosity,  and  density;  the  fluid  flow  characteristics, 
such  as  the  pressure  gradient  and  whether  the  flow  is  laminar  or 
turbulent;  and  the  geometry  of  the  system.  Hence,  of  primary  interest 
is  the  determination  or  prediction  of  the  magnitude  of  the  local  heat 
transfer  coefficient,  hcz,  or  its  dimensionless  equivalent,  the  Nusselt 
number  hcz  z/kg. 

Since  it  is  more  meaningful  to  specify  a  heat  flux,  let  us  rearrange 
equation  (15)  to  the  form 


Qgas-*-solid 

A 


=  q 


hcz  <TS  -  Tg  ) 


Substitution  of  equation  (12)  into  equation  (16)  yields 


(16) 


q 


‘Mu— 


) 


which  can  be  transformed  into  a  more  useful  form  by  equation  (11)  to 
read 


-0.0288  (Pr)1/3  (Rez)4/5  k„ 

- z - <TS  -  Tg)  . 


(17) 


The  characteristic  gas  temperature,  Tg,  is  a  troublesome  param¬ 
eter  and  deserves  some  special  consideration.  Rohseno\v  and  Choi9 
point  out  that  heat  transfer  driving  function  for  high  speed  flow  is  not 
the  free  stream  or  stagnation  temperature,  but  the  adiabatic  wall 
temperature,  Taw.  Even  though  the  velocity  at  the  wall  (surface  of 
the  plate)  is  zero,  the  adiabatic  wall  temperature,  Taw,  is  equivalent 
to  the  temperature  which  would  be  reached  if  the  wall  were  insulated. 
The  adiabatic  wall  temperature  may  be  specified  by  defining  a  recovery 
factor,  RF,  to  be 


RF  = 

where  as  previously  defined 


T*,  is  free  stream  temperature 
T0  is  stagnation  temperature. 


(18) 


^Rohsenow,  W.  M. ,  and  Choi,  H. ,  Heat,  Mass,  and  Momentum  Transfer, 
Prentice -Hall,  Englewood  Cliffs,  N,  J.,  1961. 

^Olson,  R.  M. ,  Essentials  of  Engineering  Fluid  Mechanics,  Inter¬ 
national  Textbook  Company,  Scranton,  Penna.  ,  1962. 
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Ackerman*  *  studied  theoretically  the  turbulent  boundary  layer  in  higl^ 
velocity  flow  <v*-r  a  flat  rate.  For  ordinary  high  speed  flow,  Ackerman** 
showed  that  the  fallowing  relationship  for  the  recovery  factor  is  satis¬ 
factory  RF  =  (Pr)*-  ^  .  However,  since  the  flow  conditions  in  this 
analysis  are  so  drastically  different  than  those  of  ordinary  high  speed 
flow  some  precautions  are  in  order.  The  Prandtl  number  for  the 
propellant  gates  with  which  this  analysis  is  concerned  is  2.  905; 
from  the  above  equation  RF  =  1. 43.  The  fact  that  RF  exceeds  unity 
presents  an  alarming  result.  With  RF  >1,  equation  (IS)  becomes 


aw 


T  >  T  -  T 

CO  ^  o  00 


Taw  >  TQ 

This  implies  that  the  adiabatic  wall  temperature  can  reach  values  above 
the  stagnation  temperature  which,  of  course,  is  impossible  unless 
chemical  reactions  are  present.  Therefore,  let  us  set  RF  =  1  so 
that  the  adiabatic  wall  temperature  is  equal  to  the  stagnation  tempera¬ 
ture.  Therefore,  the  characteristic  gas  temperature,  T  ,  and  the 
adiabatic  wall  temperature,  T_„„  can  be  redefined  as  the  stagnation 
temperature  Tq.  *W 


Heat  Transfer  Withi  i  the  Solid 


Subject  to  the  correct  boundary  conditions,  conduction  in  the  solid 
is  governed  by  the  Classical  Heat  Conduction  Equation  which  is 


T, 


i_  a  ts 

<  at 


(19) 


where  To  is  the  temperature  at  any  point,  x,  and  at  any  time,  t,  in  the 
solid,  ana  K  is  the  thermal  diffusivity  =  kg/p-  C_g.  Since  Figure  6 
portrays  a  cylindrical  geometry,  one's  first  thought  might  be  to  cast 
equation  (19)  into  cylindrical  coordinates.  However,  it  is  possible 
to  make  a  simplifying  assumption  enabling  one  to  solve  equation  (19) 
in  a  semi-infinite  solid  subject  to  a  boundary  condition  described  by 
equation  (17)  with  initial  solid  temperature,  T^. 


**Ackerman,  G.,  Forsch,  Gebiete  Ingenieurw,  Vol.  13,  1942. 
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Ozisik*2  gives  the  temperature  distribution  in  a  region  exterior 
to  a  cylindrical  hole  of  radius  a ,  extending  to  infinity,  and  subject 
to  a  constant  heat  flux.  From  his  analysis  it  is  possible  to  identify 
a  phenomenological  distance  e(t),  called  the  thermal  layer,  which 
represents  the  distance  beyond  which  the  initial  temperature  d  .8- 
tribution  remains  unaffected  by  the  boundary  condition.  The  esults 
of  the  analysis  show  that  the  thickness  of  the  thermal  layer  '.s  given 
by 

e  (t)  =  /5TT 

^aluminum  =  0.  20  x  10“3  inch 
f(f)brass  =  0. 18  x  10' J  inch 

provided  a  >  c  (t)  and  the  times  are  relatively  short.  Ozisik1  further 
points  out  that  for  the  same  short  times  the  above  equation  also  represents 
the  thermal  thickness  that  would  be  obtained  for  a  semi-infinite  region 
in  a  cartesian  system  with  a  constant  heat  flux  applied  at  x  =  0. 
(Admittedly,  th»*  constant  heat  flux  boundary  condition  is  a  different 
type  than  that  described  by  equation  (17).  However,  this  analysis  is 
used  only  to  demonstrate  that  the  thermal  layer  does  not  penetrate 
the  solid  to  any  significant  depth. )  Hence,  the  semi-infinite  solid 
can  be  justified  as  a  limiting  case  which  gives  results  to  the  correct 
order  of  magnitude. 

Equation  (17)  represents  a  Newton's  Law  of  Cooling  type  boundary 
condition  of  the  third  kind.  Equation  (19)  then  reduces  to 


92Ts(x,t)  1  3Ts(x,  t) 


subject  to 

q  (0,  t)  =  - 

with  an  iniKal  condition  that  Tg  (x,  0)  -  and  the  additional  boundary 
condition:  x-**«,  T§  (x,  t)-**-0.  Figure  10  portrays  a  graphical  repre¬ 
sentation  of  the  conduction  problem. 

Rohsenow  and  Choi^  give  a  solution  to  such  a  problem  and  further 
point  out  that  the  heating  of  a  thick  body  by  a  hot  fluid  at  the  surface 
approximates  this  case  during  the  early  stages  of  the  transient. 

Their  solution  is  as  follows 


^Rohsenow,  W.  M.,  and  Choi,  H. ,  Heat,  Mass,  and  Iviomentum 
Transfer,  Prentice-Hall,  Englewood  Cliffs,  N.  J. ,  1961. 

12** 

Ozisik,  M,  N.,  Boundary  Value  Problems  of  Heat  Conduction,  Inter¬ 
national  Textbook  Company,  Scranton,  Penna.,  1 9^8. 


0.0288  (Pr)  (Rez) 


(Ts  "  Tg) 
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(20) 


At  the  bore's  surface  (x=0),  the  complexity  of  equation  (20)  is  reduced 
considerably.  Although  this  equation  has  distance  into  the  solid  (x) 
and  time  (t)  as  the  independent  parameters,  the  position  along  the 
bore's  axis  (z)  is  hidden  in  hcz.  For  this  reason,  the  local  heat 
transfer  coefficient 


h 


cz 


Nuz  kg 


0.0288  (Pr)1/3  (Rez)4/5ks 


(21) 


must  be  evaluated  at  various  values  of  z.  However,  the  temperature 
atwhich  the  fluid  properties  are  evaluated  is  a  reference  temperature, 
T*,  determined  empirically  by  Eckert^  to  be 


T*  =  +  0.  50  (T0  -  Too)  +  0.  22  (Taw  -  Too)  . 

But  since  Taw  is  set  equal  to  T0,  the  above  equation  becomes 

T*  =  Tro  +  0.  50  (T0  -  Too)  +  0.  22  (T0  -  Too)  =  Tm  +  0.  72  (T0  -  Too)  (22) 


Equation  (22)  allows  determination  of  the  reference  temperature,  T*  ; 
this  value  is  then  used  to  evaluate  the  and  the  v  required  by  the 
Rez  appearing  in  equation  (21).  A  typical  calculation  for  hcz  at 
z  =  0.  030  inch  is  given  in  Appendix  D.  The  general  behavior  of  the 
heat  transfer  coefficient  as  a  function  of  position  along  the  bore's 
axis  is  given  in  Figure  11. 

Rohsenow  and  Choi^  provide  a  graphical  solution  to  equation  (20) 
in  the  form  of  a  plot  of  dimensionless  tempe rature  ,(T8  -  T^  /  Tg  -  T^)  , 
as  a  function  of  the  dimension^ ^  distance  , x/y  4 Kt  ,for  several  different 
values  of  the  dimensionless  time  »h~zy  Kf  /k  .  Since  th’S  analysis  is 
concerned  only  with  the  bore's  surface  (x  =  0  ),  Rohsenow  and  Choi's 
plot  is  reduced  to  the  values  of  the  dimensionless  temperature  along 
the  ordinate.  For  this  reason,  it  is  more  meaningful  to  plot 


^Rohsenow,  W.  M. ,  and  Choi,  H. ,  Heat,  Mass,  and  Momentum  Transfer, 
Prentice-Hall,  Englewood  Cliffs,  N.  J. ,  1961. 

13Eckert,  E.  R.  G, ,  Transactions  of  the  American  Society  of  Mechanical 
Engineers,  Vol.  78,  1956. 
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POSITION  ALONG  BORE'S  AXIS  (INCHES) 


(Ts  -  T^)/Tg  -  Tj)  as  a  function  jf  the  dimensionless  time,hcz^F  /k8 , 
as  is  done  in  Figure  ±2.  To  determine  the  solid  temperature,  Tg, 
as  a  function  of  time  at  specified  values  of  z,  the  stagnation  temperature 
and  the  heat  transfer  coefficient  (for  that  value  of  z)  are  first  obtained. 
Both  the  heat  transfer  coefficient  and  stagnation  temperature  are 
assumed  to  be  independent  of  the  solid  material.  Hence,  the  surface 
temperatures  as  a  function  of  time  for  two  materials,  aluminum  and 
brass,  are  easily  obtained  by  substituting  the  appropriate  values  of 
K  and  kg  in  the  dimensionless  parameter  h„zy/Kt/ks.  Figures  13  and 
14  are  the  bore  surface  temperature  as  a  function  of  time  for  aluminum 
alloy  707^  and  conventional  70-30  brass,  for  the  value  of  z  as  specified. 


DISCUSSION 


The  analytical  model,  developed  in  the  preceding  sections,  cul¬ 
minating  in  equation  (20),  allows  determination  of  the  bore  surface 
temperature  given  gas  flow.  To  investigate  the  performance  of  brass 
and  aluminum  alloy  cartridge  cases,  it  is  a  simple  matter  to  substitute 
the  appropriate  physical  parameters -i.  e. ,  thermal  conductivity  and 
thermal  diffusivity  -and  calculate  the  surface  temperature  as  was  done 
to  prepare  Figures  13  and  14.  Accepting  the  fact  that  aluminum  alloy 
7075  melts  in  a  range  of  890  to  1180°F  whereas  cartridge  brass  has 
a  melting  range  of  1680  «,o  1750°F,  it  is  quite  evident  that,  given  gas 
flow  for  2.0  milliseconds,  the  aluminum  melts.  The  brass  surface 
temoerature,  on  the  other  hand,  is  lower  than  its  melting  temperature. 

This  analysis  merely  indicates  that  melting  can  contribute  to 
failures  in  aluminum  cartridge  cases.  It  is  well  known  that  aluminum 
is  a  very  chemically  reactive  metal.  Hence,  there  exists  the  possibility 
c  '  reaching  a  critical  temperature  at  the  bore's  surface  which  initiates 
an  exothermic  reaction.  This  exothermic  reaction  can,  of  course, 
greatly  alter  the  previously  described  heat  flux.  There  is  also  another 
possibility.  Small  pieces  of  aluminum  may  be  torn  from  the  surface  of 
the  case  and  may  find  themselves  in  the  hignly  reactive  gaseous  plume 
exiting  the  induced  fissure.  It  is  then  possible  for  an  exothermic  re¬ 
action  to  occur  in  the  gaseous  atmosphere  which  also  alters  the  heat 
flux  to  the  surface. 
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CONCLUSIONS 


It  is  concluded  that: 

a.  Propellant  gas  flow  through  an  induced  path  (a  drilled  hole) 
in  an  aluminum  alloy  cartridge  case  does  increase  the  surface  tem¬ 
perature  in  the  path  to  the  melting  point  when  the  gas  path,  or  bore 
diameter,  is  0.0135  inch  and  the  internal  ballistic  cycle  is  normal  for 
a  5.  56  mm  cartridge. 

b.  Melting  of  the  metal  in  the  gas  path  surface  initiates  early 
(approximately  0.25  milliseconds  elapsed  time)  in  the  internal  bal¬ 
listic  cycle. 

c.  When  this  analysis  is  applied  to  a  brass  cartridge  case, 
under  the  same  conditions  mentioned  in  "a"  above  for  the  aluminum 
case,  the  melting  point  of  gas  path  surface  ie  only  approached  during 
the  entire  internal  ballistic  cycle. 


RECOMMENDATIONS 

It  is  recommended  that  the  general  analysis  to  construct  a  model 
that  describes  the  nature  of  the  aluminum  case  failure  phenomenon  be 
continued  to  include  the  effects  of: 

a.  propellant  gas  constituency 

b.  case  wall  thickness 

c.  propellant  temperature  and  pressure 

d.  size  of  the  initial  gas  path 

e.  dynamics  of  gas  path  growth 

f.  contribution  of  exothermic  chemistry  on  the  heat  flux 

g.  physical  surface  phenomenon  of  the  gas  path 

h.  additional  transport  phenomena-mass  and  momentum  transfers 

i.  influence  of  shock  fronts 

These  effects  will  be  considered  in  relation  to  the  interior  bal¬ 
listic  cycle. 
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APPENDIX  A 


System  Characteristics 


I.  Interior  Ballistics  Data 
Standard  Charge,  C  =  27.  0  grains 
Average  Gas  Pressure,  PQ  =  2,  5  x  104  psi 

II.  Propellant  Composition 


Ingredient 

Percentage 

Nitroglycerin 

9.71 

Dinitrotoluene 

0.71 

Diphenylamine 

0.90 

Dibutylphthalate 

5.96 

Nitrocellulose  (13.5%  N) 

81.96 

Moisture  and  volatiles 

0.85 

Residual  solvent 

0.29 

Calcium  carbonate 

0.46 

Sodium  sulfate 

0.07  (nominal) 

III.  Propellant  Product  Gas  Properties 

Adiabatic  flame  temperature,  TQ(*4)  =  2800 °K=  5040°  R 

Gas  constant,  R  =  64,372  ft-lbf/lb“R 

i  m 

Moles  of  gas  produced,  n^4)  =  0.  0423  gram-moles/gram  of  propellant 

Specific  heat  at  constant  volume,  Cv(*^  =  0.344  BTU/lbm°R 

Specific  heat  at  constant  pressure,  Cp  =  0.427‘"  BTU/lbm°R 

Ratio  of  specific  heats,  y(^4)  =  Cp/Cv  =  1.24 

Molecular  weight,  MW  (16)  =  23.62 

Absolute  viscosity,  =  4.72  x  10“5  lbm/ft-sec 


l4Unpublished  Data,  "Smokeless  Powder  for  Small  Arms  Ammunition, 
Nominal  Characteristics  of  Standard  Products,  "  Olin  Corporation, 
New  Haven,  Conn. ,  Oct.  1966. 

l^Price,  E.  W. ,  et  al. ,  "Theory  of  Ignition  of  Solid  Propellants,  " 
AIAA  Journal,  Vol.  4,  July  1966. 

^Stiefel,  L. ,  and  Hody,  G.  L.,  "The  Composition  of  the  Exhaust 
Products  of  Military  Weapons,  "  Technical  Report  R-1948,  March 
1970,  Frankford  Arsenal,  Philadelphia,  Penna. 
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Thermal  conductivity,  k  = 

O 

Prandtl  number,  Pr  =  2. 905 
Recovery  factor,  RF  =  1.43 
Propellant  gas  composition  16) 

Compound 

H  2 
H20 
CO 
N2 

co2 


2.  5  x  10-2  BTU/hr-ft  °R 


Composition  in  Mole  Fraction 

0. 1531 
0. 1968 
0.4419 
0. 1078 
0. 0986 


*^Price,  E.  W.,  et  al. ,  "Theory  of  Ignition  of  Solid  Propellants,  11 
AIAA  Journal,  Vol.  4,  July  1966. 

^Stiefel,  L. ,  and  Hody,  G.  L.,  "The  Composition  of  the  Exhaust 
Products  of  Military  Weapons,  "  Technical  Report  R-1948,  March 
1970,  Frankford  Arsenal,  Philadelphia,  Penna. 
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APPENDIX  B 


Determination  of  Intermediate  Mach  Number 


Bore  length,  L 
Bore  diameter,  D 
Friction  factor,  f 
Inlet  Mach  number,  M 


=  0.  062  inch 
=  0. 0135  inch 
=  0.002 
=  0.  865 


1.  The  bore  length  required  for  the  Mach  number  to  change  from 
0.865  to  0.90. 


=  .0367  -  .0174  =  .0193 


L 


(.  013")  [.  jL'W 
(4)  {.  002) 


0.  032  inch 


The  Mach  number  increases  to  0.  90  after  0.  032  inch. 

2.  The  bore  length  required  for  the  Mach  number  to  change  from 
0.90  to  0.95. 

4fL  _  /  4fL<\  /  4flA  AflA  / _4fJL.\ 

D  \  °  /M1  "  \  °  /M2  \  °  A  900  \  D  /  .  950 

=  .0174  -  .00395  =  .01345 


(.  0135)  (.01345) 
(4)  (.002) 


0. 0226  inch 


The  Mach  number  increases  from  0.  865  to  0.  95  after  0.  032  inch 
+  0,  0226  inch  =  0.0546  inch. 


3.  The  bore  length  required  for  the  Mach  number  to  change  from  0.  95 
to  1. 00. 
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4fL  _  /  4fL \  / 4fL.\  _  /  4flA  f 4fL  | 

D  "  \%1  “  \D/M2"  V  D  /. 95  "  VD/l. 


00 


=  .00395  -  0  =  .00395 


t  .  (.0135)  (.00395)  _  „ 

L  =  " (4)  (.0  0?) -  -  °*0066inch 


The  Mach  number  increases  from  0.865  to  1.00  after  0.032  inch 
+  0. 0226  inch  +  0.  0066  inch  =  0.0612  inch. 


Summary  for  f  =  0.002 


z 

M 

0 

0.865 

0.032 

0.90 

0.0546 

0.95 

0.0612 

1.00 
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APPENDIX  C 


Calculation  of  Gas  Properties  at  Intermediate  Locations  (f=0.002) 


Chamber  conditions 

T0  =  5040 °R 
PQ  =  2.5  x  104  psi 
po  =  11.09  lbm/ft5 

1.  z  =  0  (inlet  conditions) 

Isentropic  Process 


M 

=  .865 

T 

00 

=  4624 °R 

(equation  (4)) 

Pm 

00 

=  16025  psi 

(equation  (5)) 

P  X 

=  7.  75  lbm/ft3 

(equation  (6)) 

vro 

-  2981  ft/sec 

(equation  (?)) 

Po 

=  2.5  x  104  psi 

(given) 

To 

:  5040 °R 

(given) 

2,  z  =  0.062  (exit  conditions) 

Adiabatic  Process 


M 

To 

Too 

P_ 

P* 

P* 

Poo 

Px 


1.0 

5040  °R.  (given) 

4500  °R  (equation  (4),  M  =  1) 


1.  173 

16025 
1.  173 


=  13661  psi 


(64.  372)  (4500) 
6.79  lbm/ft3 


(gas  tables) 
(equation  (3)) 


P0  =  (1.  795)  (13661)  (equation  (5),  M  =  1) 

PQ  =  24521  psi  _ 

v  =  1  y/Jl.24)  (2072.  7)  (4500)  (equation  (7),  M  =  1) 

v  =  3400  ft/sec 
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M 

— 

0.95 

T 

xo 

= 

5040 °R 

(given) 

T 

00 

= 

4547 °R 

(equation  (4),  M  =  .  95) 

P 

P* 

= 

1.0583 

(gas  tables) 

P 

= 

(1.  0583)  (13661) 

P 

CO 

= 

14458  psi 

Poo 

= 

(14458)  (144) 
(64.372)  (4547) 

(equation  (3)) 

Poo 

= 

7.  11  lbm/ft3 

V  =  .95/(1. 24)  (2072,  7)  (4547)  (equation  (7),  M  =  .  95) 
v  =  3247  ft/sec 

PQ  =  (1.  705)  (14458)  (equation  (5),  M  =  .  95) 

Pq  =  24650  psi 
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